A tunable source of correlated atom beams 
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We use a one-dimensional optical lattice to modify the dispersion relation of atomic matter waves. 
Four-wave mixing in this situation produces atom pairs in two well defined beams. We show that 
these beams present a narrow momentum correlation, that their momenta are precisely tunable, and 
that this pair source can be operated both in the regime of low mode occupancy and of high mode 
occupancy. 
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In quantum optics, existence of mechanisms to pro- 
duce photon pairs, such as parametric down conversion, 
enabled the realisation of several fundamental experi- 
ments on quantum mechanics. For example, the vio- 
lation of Bell's inequalities |T] or the Hong-Ou-Mandel 
effect [2j reveal the surprising properties of quantum cor- 
relations in entangled photon pairs. These fascinating 
properties have found applications in quantum informa- 
tion and communications [3j. In analogy to photon pairs, 
there have been several recent demonstrations of corre- 
lated atom pairs production [4]-[9]. In particular, momen- 
tum correlations of spatially separated samples is an im- 
portant requirement for the demonstration of an atomic 
EPR state plOl [TT] and the violation of Bell's inequal- 
ities. Such momentum correlations were demonstrated 
for atom pairs produced by molecule dissociation [4J or 
by spontaneous four-wave mixing in free space through 
the collision of two Bose-Einstein condensates (BEC's) 
[5l [12] . In these experiments the pairs which were pro- 
duced lay on a spherical shell. This geometry is disad- 
vantageous because many spatial modes are populated, 
and if one wishes to use Bragg diffraction to manipulate 
and recombine the pairs on a beam splitter [TOl [13] , the 
vast majority of the pairs is unuseable. 

On the other hand, if pair production can be concen- 
trated in a small number of modes, experimenters can 
make more efficient use of the generated pairs. One can 
then chose to work either with low mode occupation, the 
well separated pair regime, or with high mode occupa- 
tion, refered to as the squeezing regime in Ref. [M]. An 
example of twin beams generated in the latter regime is 
described in Ref. [7J . The squeezing regime is well suited 
to the study of highly entangled multiparticle systems, 
and for investigations of atom interferometry below the 
standard quantum limit [151 HSj- The pair source we 
study in this letter can be operated in both regimes. We 
use atomic four- wave mixing in a one-dimensional (ID) 
optical lattice, which results in production of atom pairs 
in two well defined beams, as proposed in Ref. [17J and 
demonstrated in Ref. [18j. We show that these beams 
present a narrow momentum correlation, that their mo- 
menta are precisely tunable, and that we can control their 
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FIG. 1. (Color online) (a) ID pair creation process in an 
optical lattice with period Xiatuce / The dispersion relation 
in the first Bloch band (green solid curve) allows scattering 
of atoms from a BEC with quasi-momentum /co (open red 
circle) in the lattice frame into pairs with quasi-momenta ki 
(filled orange circle) and /c2 (filled blue circle), so that phase- 
matching conditions given by energy and momentum conser- 
vations are fulfilled. The example here is for a lattice depth 

Vb = 0.725 Erec and /Co = -0.65 /Crec, with k^ec = ^^TI / XlatUce 

the recoil momentum and ^rec = h^k^ec/'^'n^ = /i x 44 kHz the 
recoil energy. 

(b) Vertical single-shot momentum distribution (integrated 
over the total transverse distribution) measured for these con- 
ditions. The three main peaks correspond to the initial BEC 
and to the macroscopically populated beams centered at ki 
and /c2, which are mainly projected in the first Brillouin zone 
(in white) when the lattice is switched off. As expected, small 
diffraction peaks at ko + 2/crec and /c2 — 2/crec are also visible, 
due to the proximity of ko and /c2 to the band edge. 



intensities. 

In atom optics, four- wave mixing corresponds to scat- 
tering into new momentum classes subject to energy and 
momentum conservation. In a wave picture, the conser- 
vation requirements can be thought of as phase matching 
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FIG. 2. (Color online) Experimental setup and sequence: 

(1) Initially, a BEC of metastable helium is trapped in a ver- 
tical optical potential with a 43 pm waist. 

(2) An optical lattice is suddenly applied in the presence of 
the trap. It is tilted by 7° with respect to the trap axis, and 
is focused on the BEC with a 200 pm waist. 

(3) After dipole trap and optical lattice switch off, the cloud 
expands and falls on the 3D resolved single atom detector. 
Given the values of the vertical and transverse Thomas- Fermi 
radii (0.5 mm and 3 pm), the arrival time and position reflect 
the 3D momentum distribution, provided the momenta are 
well above 3 x 10~^ krec along z and 2 x 10~^ krec transversely. 



conditions. The presence of an optical lattice modifies the 
free-space atomic dispersion relation and therefore, for a 
range of initial quasi-momenta /cq [19j, the ID scattering 
event 2 /cq ^ /ci + /c2 is allowed, as shown in Fig. [ija). 
Thus, beginning from a BEC at /cq, atom pairs are spon- 
taneously generated along the lattice axis with well de- 
fined quasi-momenta ki and /c2. We will refer to this 
process as four-wave mixing, but it can also be thought 
of as a special case of a dynamical instability [2Q[ [21], 
which was studied in the context of coherence [2J, ^ 
and atomic [24] losses appearing for a BEC moving in a 
lattice. 

The experiment is performed on ^He atoms in the 
rrix = 1 sublevel of the 2^5'i metastable state. The ex- 
perimental setup and sequence are shown in Fig. [2] After 
evaporative cooling in an elongated, vertical dipole trap 
with frequencies = 1.5 kHz and Uz = 6.5 Hz [25j, we 
produce a BEC (or more precisely a quasi-BEC [26j) with 
about 10^ atoms. We then apply an optical lattice with a 
depth Vq = 0.725 ^rec- This lattice is tuned 19 nm to the 
blue of the 1083 nm 2^6'i - 2^P transition of helium. It 
is formed by two counterpropagating 17 mW beams with 
200 fim waists and whose relative detuning 5iy can be var- 
ied using acousto-optic modulators. We thus control the 
value of ko/krec = hdv/AE^Q^^ the BEC's momentum in 
the lattice frame. The lattice is held on for a duration 



Tf, = 2 ms, and suddenly switched off, simultaneously 
with the optical trap. In order to avoid perturbation of 
the cloud during time-of-flight due to residual magnetic 
gradients, we apply an RF pulse that transfers 50% of the 
atoms in the field insensitive rrix = sublevel ^25]. The 
atoms remaining in rrix = 1 are subsequently removed by 
a strong magnetic gradient. After a 307 ms mean time 
of flight, the rrix = atoms fall on a micro-channel plate 
detector, which permits 3D reconstruction of the atomic 
cloud [27j. 

As shown in Fig. [ijb), we detect three main density 
peaks after the time of flight. The higher one is the 
initial BEC. The two other peaks are formed by the 
atoms scattered into the momentum classes centered in 
ki and /c2, whose values are consistent with those ex- 
pected from the phase-matching conditions illustrated in 
Fig. [ija). Since the periodic potential is switched off 
abruptly, the Bloch states of wavevectors /cq, ki and k2 
are projected onto plane waves, mainly in the first Bril- 
louin zone due to the low lattice depth. Each of the 
macroscopic beams at ki and /c2 contains on the order 
of 10^ detected atoms, which we estimate to correspond 
to about 2 X 10^ atoms per beam. We also detect some 
atoms between the beams, which result from scattering 
into excited transverse modes [28j. Due to the low over- 
lap between the transversely excited states and the initial 
wave function, this transverse excitation is far less effi- 
cient than the previously described ID process. In addi- 
tion, scattered atoms can undergo secondary scattering 
also contributing to the background between the density 
peaks. 

In the following, we focus on the two beams. Using 
them for quantum atom optics experiments or for in- 
terferometry will require recombining them. It is there- 
fore crucial to know the width of their correlation. From 
the momentum distribution n(k), we computed the nor- 
malised second-order cross correlation function, 



{n(k)n(k')> 
(n(k)) (n(k')) 



(1) 



where k belongs to beam 1 and to beam 2. The BEC is 
not exactly at rest in the optical trap, but exhibits shot 
to shot momentum fluctuations of order 10~^ /crec- We 
correct for these fluctuations by recentering separately 
the single shot momentum distributions n(k) around ki 
and /c2, using the shift obtained from Gaussian fits to the 
peak at ki and to the diffraction peak at k^ + 2/crec- As 
shown in Fig. [3]^a,b), this correlation function exhibits a 
peak for kz — ki and k'^ k2' The presence of this peak 
indicates that the two atomic beams are indeed corre- 
lated. 

We wish to determine the number of modes present in 
each beam, and how many of these modes are correlated. 
We therefore examine the local second-order correlation 
function of a single beam, (k, k'), which is obtained 
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FIG. 3. (a,b) Cuts along {y,z) of the integrated, normal- 
ized cross correlation function of the two beams, ^^'*(Ak) — 
f dki ^^'*(ki, kj + Ak). The integration over the momen- 
tum distribution ki is performed on a box with dimensions 
Lk^ — Lky = 0.4 krec and Lk^ = 5 x 10~^ krec centered on 
beam 1, ki + kj = (/ci + fe)ez, and the cuts have a thick- 
ness 10~^ krec (1-5 X 10""^ krec) aloug z {x and y). The 
bunching, due to the correlation between the cross beams, 
has a longitudinal (transverse) width ac,z = 1.8 x 10~^ krec 
(o'cy = 1.6 X 10""*^ krec)- (c,d) Cuts aloug (y^z) of the in- 
tegrated, normalized local correlation function of beam 1, 
g^^\Ak) = / dki^[^^(ki,ki + Ak). The integration zone 
is the same as for the cross correlation, and the cuts have a 
thickness 2.5 x 10~^ krec (0.1 krec) along z (x and y). The 
bunching, due to HBT effect, has a longitudinal (transverse) 



width ai^z = 3.7 x 10 krec {(yi,y = 1-3 x 10 fc, 



Cuts 



along X (not shown here) have same widths and amplitudes 
as cuts along y. These correlation functions are calculated 
using 850 experimental realisations, with ko = —0.65 krec and 
a lattice depth Vo = 0.725 ^rec- 



as in Eq. [T]but with both k and k' belonging to beam 1. 
This correlation function, plotted in Fig.|3jc,d), exhibits 
bunching for k'^ — kz — /^i, due to density fluctuations 
(as in the Hanbury Brown and Twiss effect [22] )• Sim- 
ilar bunching is observed at /c2. If we suppose that the 
widths of the local correlation define the size of a sin- 
gle mode, we can compare them to those of the den- 
sity (longitudinal rms: 4 x 10~^ /crec, transverse rms: 
4 x 10~^ Kec)' We see that about 10 longitudinal and 3 
transverse modes are populated. Thus the mode popu- 
lation is, very roughly, 70 atoms/mode. For comparison, 
in the case of free space four- wave mixing [30], starting 
from a similar initial BEG, 10^ modes were populated, 
with only about 0.02 atom/mode. 
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FIG. 4. Measured mean momenta ki and k2 of the twin beams 
(black dots, in units of krec) as a function of the initial BEG 
momentum ko (in the lattice frame). The solid line shows the 
phase-matching curve expected in the absence of interactions, 
while the dashed line includes the mean field (see Eq.|2|. 



It appears in Fig. [3] that, while in the transverse direc- 
tion, the cross and local correlations have similar widths, 
the cross correlation is 5 times broader than the local one 
along the vertical axis: each mode is correlated with sev- 
eral modes of the other beam. If one uses two such beams 
as inputs to a beam splitter, this broadening amounts to a 
loss of coherence, and the interference contrast would be 
reduced. We emphasize that the observed widths may be 
broadened by other effects, and so their numerical ratio is 
not exactly equal to the number of correlated modes. For 
the cross correlation, the observed width is broadened by 
the fact that the vertical source size is not entirely negli- 
gible [31j. For the local correlation, we estimate that the 
finite vertical resolution of the microchannel plate detec- 
tor contributes notably to the observed width. This reso- 
lution comes about because the surface which defines the 
atom arrival time is not flat but consists of tilted channels 
which intercept the atoms at different heights. The width 
shown in Fig. (sjd) is consistent with this interpretation. 

The use of an optical lattice permits control over the 
output beam momenta. Changing the detuning Su be- 
tween the lattice beams results in varing the value of ko. 
In Fig.[4j we plot the mean vertical momenta ki and k2 of 
both beams, measured for different ko, as well as the ex- 
pectation (continuous line) based on the phase matching 
conditions illustrated in Fig.jlja). We obtain a fair agree- 
ment over a large range, even though the continuous line 
presents a small shift in comparison to the data points 
and does not reproduce the observed shape for high val- 
ues of ko. However, as already observed for four- wave 
mixing in free space [32j, phase matching conditions can 
be influenced by mean field effects. A simple correction 
to the phase matching curve is found just by adding the 
mean field to the energy conservation condition: Since 
the two atoms of a scattered pair are distinguishable from 
the atoms of the initial BEG, the mean-field energy expe- 
rienced by each of them is not gno (with g = 47r fi^a/m, 
a and m the scattering length and the mass of He* and 
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FIG. 5. (Color online) Dependence of the population of beam 
1 (orange filled circles) and beam 2 (blue open circles) on 
the lattice duration Tl for ko = —0.67 krec and for a lattice 
depth Vb = 1-05 ^rec (This depth is a bit higher than for 
the previously described data, which increases the process 
efficiency). The gray line is an exponential fit of the detected 
population in beam 2 for Tl < 0.3 ms, which gives a time 
constant of 0.103 ms and an offset of 11.5 detected atoms. 
This offset is due to the small thermal part of the source 
cloud with quasimomenta ki and A;2. Inset: same data with 
linear scale. 



no — 10^^ atoms/cm^ the EEC's density), but 2^no, so 
that the energy conservation condition reads: 

2E{ko) + 2gno = E{ki) + E{k2) + 4^no (2) 

where the energy E{k) is given by the dispersion relation 
in the first Bloch band of the lattice without interaction. 
As seen in Fig. [4] (dashed line), this correction leads to 
very good agreement with the experimental data, and 
accounts for the shift of the phase-matching curve and 
the change of its shape. A more exact calculation of the 
phase-matching conditions, inspired by [20j, confirms the 
accuracy of Eq.[2]in our experimental conditions and will 
be given in [33j. 

Another degree of freedom results from the fact that 
pair creation only takes place while the lattice is on. We 
can thus tune the beam populations with the lattice du- 
ration Tl- In the example of Fig. [5] these populations 
increase exponentially with T^ during a few hundreds of 
ps, and then reach a plateau. This saturation could be 
explained by several mechanisms such as the decrease of 
spatial overlap between condensate and scattered beams 
[18j, multimode effects [34J and secondary scatterings 
from the beams. Condensate depletion is at most about 
20 %, and should be of little importance in the satura- 
tion. For small T/,, there is no discernable population 
difference between both beams. By contrast, we observe 
that at large Tl the population of beam 1 is almost twice 
that of beam 2, a phenomenon also noticed in ^8j. This 
may be due to /c2 being in a dynamically unstable region 
while atoms with quasi-momentum ki can only undergo 
secondary scattering to excited transverse modes. 



Figure [5] shows that as long as the pair production rate 
is small (lattice durations less than 0.5 ms), losses due 
to secondary scattering processes are negligible and the 
resulting beams should contain strongly correlated pairs. 
In an attempt to verify a non-classical correlation, we 
examined number difference between the two beams. By 
selecting a narrow transverse window in the two peaks, 
we do indeed observe a sub-Poissonian difference [7, 30j, 
with a variance about 10% below the shot noise. The size 
of this window is nearly an order of magnitude narrower 
than that of the transverse correlation function, a result 
that is unexpected and which we plan to investigate in 
future experiments. 

To conclude, we have demonstrated an efficient pro- 
cess for the production of correlated atom pairs. We 
have some degree of control over both the final momenta 
and the intensity of the correlated beams, and we have 
characterized the width of the correlation in momentum. 
Thus we think that this source is a good candidate for 
use in multiple particle interference experiments both in 
the regime of well isolated pairs [llj and in the regime of 
large occupation numbers [TO] . 
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